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♦Present  address: 


Abstract 


A  theoretical  investigation  of  negative-ion  formation  from  positive- 
ion-surface  scattering  is  presented  from  a  unified  point  of  view.  Based 
on  the  time-dependent  Anderson-Newns  model,  the  correlation  energy  U  is 
seen  to  play  an  important  role  in  the  two-electron  transfer  process. 
Calculations  of  the  probability  of  negative-ion  formation  are  in  good 
agreement  with  experiments  on  the  conversion  of  H+(D+)  to  H"(D")  by 
scattering  from  a  cesiated  W(100)  surface. 


Charge-exchange  processes  arising  from  monoenergetic  ion  beams  scattered 

1-5 

from  solid  surfaces  have  recently  been  the  subject  of  much  experimental 

and  theoretical^"^  interest.  Most  of  these  studies  focus  on  the  process 

of  ion  neutralization,  involving  the  transfer  of  a  single  electron  from  the 

surface  to  the  singly-charged  ion. -3,6-10  However,  owing  to  the  increase 

of  applications  involving  the  transfer  of  two  electrons  in  a  variety  of 
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situations,  such  as  plasma  fusion  *  and  the  generation  of  high-energy 
neutral  beams, **  attention  has  also  been  directed  to  processes  of  negative- 
ion  formation.^’ 

For  ion  neutralization,  the  majority  of  theoretical  work  has  been  based 

14 

on  the  Anderson-Newns  model,  originally  proposed  for  the  explanation  of  local¬ 
ized  magnetization  in  transition  metal  alloys  and  subsequently  applied  to  various 
other  problems  such  as  chemisorption  on  metals^  and  mixed  valence  in  solids. 

In  these  applications  of  this  model,  the  correlation  energy  U,  arising  from  the 
Coulomb  repulsion  between  the  two  electrons  of  opposite  spin  in  the  same  discrete 
level,  plays  a  crucial  role.  However,  in  studies  of  ion  netrualization  (or  the 
reverse  process  of  atom  ionization),  this  important  quantity  is  either  implicitly 
assumed  to  be  Infinite  or  completely  ignored.  ~  Such  approaches  are  thus 
incapable  of  accounting  for  negative  ion  formation:  U  *  «  completely  supresses 
the  transfer  of  a  second  electron,  while  if  only  single-electron  transfer  is 
considered,  U  Is  Irrelevant. 

13 

Previous  theoretical  work  on  negative-ion  formation  has  centered  on 
the  time-dependent  width  and  shift  of  the  valence  level  of  the  Impact  Ion, 
where  the  time  dependence  Is  due  to  the  motion  of  the  Ion.  It  was  assumed 
that  first  neutralization  occurs  via  a  nonresonant  Auger  process,  while  the 


nuclear  motion  of  the  ion  can  subsequently  bring  the  valence  level  into 
resonance  with  a  band  state  and  thus  effect  the  transfer  of  the  second 
electron.  The  quantity  U  again  does  not  play  an  important  role  in  this 
theory. 

In  the  present  work,  we  use  the  Anderson-Newns  model  to  examine  the 

effects  of  U  on  the  two-electron  transfer  process  in  relation  to  the  dynamics 

of  the  nuclear  motion  of  the  impact  ion.  This  represents  the  first  time 

that  charge-exchange  processes  in  surface  scattering  have  been  considered 

from  a  unified  point  of  view.  We  shall  see  that  the  negative-ion  formation 

depends  crucially  on  the  finite  value  of  U,  and  in  fact,  recent  experimental 
5 

results  can  be  explained  in  terms  of  our  theory. 

The  time-dependent  Anderson-Newns  Hamiltonian  including  the  correlation 
energy  term  is  given  as 

H(t)  =  Hq  +  V(t), 

where 

H0  ■  I'dcdoCdo  *  l  ‘kCka  Cko  *  U"dc"d,-c 

V(t)  =  kI[VM(t)c-oCko  ♦  V*ki(t)c^]. 

The  Indlcies  d  and  k  denote  the  valence  state  of  the  Impact  ion  and  the 
conduction  band  states  of  the  solid,  respectively,  o  is  the  spin  Index  and 
ndo  *  CdoCdc*  The  ^nteract*on  H#n>iltonian,  V(t),  responsible  for  the  electron 
transfer  between  the  band  and  valence  states.  Is  the  only  explicitly  time- 
dependent  part.  The  motion  of  the  projectile  ion  can  be  phenomenologically 
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taken  into  account  by  using  the  specific  form 
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=  ‘'kd6' 


i 


where  x,  the  sole  dynamical  parameter  in  our  model,  is  directly  proportional 
to  the  normal  impact  velocity.  Hence  x  controls  the  duration  of  the  bound- 
continuum  interaction,  while  (time-independent)  determines  its  strength. 

To  lowest  order,  the  perturbative  solution  for  the  time-evolution 
operator  T  in  the  interaction  picture,  which  contributes  to  the  S-matrix 
for  negative-ion  formation,  is  given  as 

T(2)(t,-»)  =  (-i)2  /  dt'V(t’)/  dtMV(t"), 


(5) 


where 

V(t)  ■  elHQtV(t)  e^V  (6) 

is  V(t)  expressed  in  the  interaction  picture.  At  t  *  -«,  the  ion  is  taken 
to  be  infinitely  far  from  the  surface,  and  t  *  0  is  the  instant  of  impact. 

It  is  apparent  from  Eq.  (3)  that,  since  negative-ion  formation  involves 
ultimately  the  transfer  of  two  electrons,  only  terms  of  even  powers  in  V 
in  the  expansion  of  T  contribute  to  the  S-matrix. 

Figure  1  illustrates  schematically  the  electron-transfer  process  to 
lowest,  i.e.,  second,  order  as  contained  In  Eq.  (S).  Physically,  every 
factor  of  V(t )  corresponds  to  the  transfer  of  one  electron  from  the  band 
level  to  the  level  of  the  Ion,  or  the  conjugate  process  of  electron 
transfer  from  to  c^.  Therefore,  the  contribution  to  the  S-matrix  from 


4 


higher-order  terms  of  T  can  be  diagrammatically  generated  in  a  straight¬ 
forward  manner.  Such  terms  (neglected  in  this  work),  however,  may  involve 
those  final  states  in  which  the  solid  becomes  electronically  excited,  with 
levels  higher  than  the  Fermi  level  occupied.  In  general,  low-ordered  pro¬ 
cesses  are  favored  by  high  impact  velocities  (large  x). 

To  lowest  order,  the  time-dependent  probability  for  negative -ion 
formation  is  then  given  as 

eF  eF 

P(t)  «  /  dc  /  d£,p(e)o(c,)|<k^:iT(Z)(t,-)|I>|2,  (7) 

eL  CL 

where  and  are  the  Fermi  energy  and  conduction  band  edge  of  the  solid, 
respectively,  p(e)  is  the  density  of  states  of  the  band,  and  the  initial 

and  final  states,  |I>  and  { k^k * >  respectively,  are  described  in  Fig.  1. 

The  matrix  element  in  Eq.  (7)  can  be  evaluated  by  using  Eqs.  (2)-(6)  to  give 

the  result 


..  .;(2)/.  \ .  _  l r  1  2x+i (U-e-e 1 )]t  f  1  A  1  1 

<k+k;  |T  '(t,-»)|I>  -  (x^T-lp-J ’ 


t  s  0 


(8a) 


V2  1 

T  X+HU-e-e1  JTT 


(x-ie  +  X-ie1)' 


V2  1 

T  xTiru“c-cH/2 


+  2V2X^ 


(x^Hx-Ku-e1)] 


(> 


e-Cx-KU-e'JH 

-  fe-[X-1(U-c)t  J 

Se,Z)[x-1(U-c)]  i 


t  >  o 


(8b) 
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where  e  =  ek  -  ed  and  c 1  =  ek,  -  ed  with  ed  set  as  the  zero  of  energy,  and 

1/  =  v,  .  is  assumed  to  be  independent  of  energy.  As  a  first  approximation,  we 
kd 

have  also  assumed  that  U  and  ed  are  constant  within  the  collision  region.  Thi 
latter  assumption  is  not  expected  to  affect  the  results  qualitatively. 

We  now  turn  to  calculations  of  the  probability  P(t)  of  negative-ion 
formation  and  its  limiting  value  at  t  -*•  «*,  P(«),  for  various  choices  of 
the  parameters  U  and  >..  P(«)  represents  the  experimental  observation  of 

outgoing  negative  ions  after  the  scattering  event  is  complete,  and  P(t) 
reflects  the  behavior  of  the  transient  states.  The  following  fixed 

2 

numerical  values  were  adopted:  =  -10  eV,  =  -0.06  eV  and  l  =  zV  = 

0.43  eV,  where  we  have  assumed  the  band  to  have  no  important  structure 
so  that  the  density  of  states  p  can  be  taken  as  constant  in  the  integra¬ 
tion  over  the  band  in  Eq.  (7).  Although  the  last  two  numbers  have  been 

g 

used  previously  in  a  study  of  charge  transfer  in  the  Na/W(100)  system, 
they  have  no  special  significance  in  our  present  work,  which  is  to 
investigate  j_n  general  the  effects  of  the  variation  of  x  and  U,  especially 
the  latter. 

Intuitively,  we  expect  the  probability  P(«)  to  be  small  when  the 
repulsive  correlation  energy  U  is  large.  However,  P(«)  also  depends  on 
x,  i.e.,  on  the  velocity  of  the  impact  ion.  As  shown  In  Fig.  2,  there  is 
a  peak  for  P(°>)  at  a  small  value  of  x  for  each  U.  The  explanation  for  this 
is  that  since  small  x  implies  a  long  duration  of  interaction  whereby  from 
the  uncertainty  principle  the  resonance  requirement  is  stringent,  it  is 
impossible  for  the  second  electron  In  the  solid  to  overcome  the  barrier  u. 

On  the  other  hand,  large  x  can  ease  the  resonance  requirement  —  energy 
conservation  can  be  violated  for  short-duration  processes  —  but  It  also 
limits  the  actual  time  available  for  electron  transfer,  resulting  In  a 
small  P(»).  As  a  consequence,  for  each  U  there  exists  an  optimal  value  of 


;■ ,  ?.  ,  for  which  the  probability  attains  a  maximum.  Moreover,  > 
m  m 

increases  as  U  is  increased,  due  to  the  fact  that  a  shorter  inter¬ 
action  time  is  favorable  for  the  second  electron  to  transfer  non- 
resonantly  as  U  becomes  larger. 

The  close  relationship  between  resonant  electron  transfer  and  the 
ion  velocity  becomes  obvious  when  we  look  at  the  time  evolution  of  the 
probability  of  negative-ion  formation.  Figure  3  displays  P(t)  for  various 
values  of  U  with  fixed  x.  Each  probability  curve  has  a  peak  at  very  short 
time  tm  ('  0.2  femtoseconds).  We  see  that  tm  is  smaller  for  larger  U,  in 
accordance  with  the  arguments  given  above.  In  addition,  the  degree  of 
transient  negative-ion  formation,  measured  by  the  ratio  P(tm)/P(»),  is  more 
pronounced  for  larger  U.  Figure  4  provides  yet  another  manifestation  of 
the  striking  transient  behavior;  namely,  although  there  is  an  optimal  xm 
at  each  U  for  the  limiting  value  of  the  probability,  P(«),  this  is  not 
necessarily  the  case  in  the  transient  region. 

To  test  our  theory,  we  shall  compare  it  with  measurements  on  the 
conversion  of  H  (D+)  to  H  (D  )  by  scattering  from  a  cesiated  W(100) 
surface  at  different  gra2ing  angles  e.5  For  this  purpose,  we  identify 
x  as  vcoss,  where  v  is  the  magnitude  of  the  velocity,  and  phenomeno¬ 
logically  introduce  a  velocity-dependent  interaction,  a  *  0.43  exp(-O.Olv) 
(In  the  unit  of  eV),  to  account  for  the  loss  of  particles  due  to  penetra¬ 
tion  Into  the  surface.  The  variation  of  e  is  thus  equivalent  to  the 
variation  of  the  normal  Impact  velocity  A.  Our  results,  given  In 
Fig.  5,  are  In  qualitative  agreement  vlth  the  experimental  ones,5 
where  for  all  Incident  enerc  «  th'  .unversion  probability  goes 
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through  a  maximum.  Quantitative  comparisons  have  not  been  attempted  since 

precise  information  on  critical  parameters,  especially  t,  is  still  lacking. 

For  the  cesiated  W(100)  surface,  among  other  complications  leading  to 

unreliable  data  for  parameters  is  the  theoretical  evidence  of  a  lowering 

18  19 

of  the  work  function  by  multiple  dipole  formation.  ’ 

In  this  work  we  have  demonstrated,  through  varying  the  dynamical 

conditions  of  the  impact  ion,  the  significance  of  the  correlation  energy  U 

in  negative-ion  formation  from  positive-ion-surface  scattering.  Though  U 

in  general  decreases  the  probability  for  n^ative-ion  formation,  one  can 

always  exploit  the  experimentally  controllable  dynamical  conditions  (varying 

v  and  e)  to  achieve  an  optimal  result  for  a  given  system.  Moreover,  there 

may  even  be  the  possibility  of  exploiting  the  characteristic  transient  behavior 

of  P(t),  since  for  finite  U,  P(tm)  is  always  larger  than  (»)  except  for  very 

large  x.  For  very  small  x,  on  the  other  hand,  our  perturbation  approach  may 

not  yield  correct  results,  since  the  long  interaction  times  then  allowed  may 

require  higher-order  processes  than  the  second-order  one  be  considered.  Our 

results  have  been  shown  to  be  in  good  qualitative  agreement  with  experiments. 

A  more  elaborate  calculation  is  needed  which  takes  into  account  the  lowering  of  the 

18  19 

valence  level  of  the  ion  near  the  surface  ’  is  needed  for  quantitative 
comparision  with  experiments. 
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Figure  Captions 

Fig.  1.  Schematic  diagram  of  electron  transfer  in  negative-ion  formation. 
Gy  is  the  upper  and  eL  the  lower  edge  of  the  band,  (a)  Initial 
state  |I>:  valence  state  empty,  band  filled  up  to  the  Fermi 
level  cp.  (b)  Intermediate  states  |k+>  and  |k|>:  states  cor¬ 
responding  to  the  neutralized  atom;  one  electron  transferred 
from  the  ck  or  e^,  level  to  the  level.  The  arrows  denote 
the  spin  states  of  the  electrons,  and  the  solid  and  hollow 
circles  represent  electrons  and  holes,  respectively,  (c)  Final 
states  1 k  k ' > :  negative-ion  states;  two  electrons  transferred 

i  > 

to  the  level. 

Fig.  2.  P(»)  vs  x  for  various  values  of  U.  Energy  is  in  the  unit  of  eV. 

Fig.  3.  P(t)  vs  t  for  various  values  of  U  with  fixed  x.  As  U  increases, 
the  characteristic  short-time  behavior  becomes  more  pronounced. 
Energy  is  in  the  unit  of  eV,  and  time  is  in  the  unit  of  6.59  x 

l(f16  sec  (J^) 

Fig.  4.  P(t)  vs  t  for  various  values  of  x  with  fixed  U.  The  units  are  the 
same  as  in  Fig.  3. 

Fig.  5.  P(«)  vs  e,  the  incident  angle  of  impact.  vi<v2<v3<v4* 
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